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This  report  describes  an  inhouse  effort  conducted  in  the  Mechanics 
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(MB),  Air  Force  Materials  Laboratory,  Air  Force  Wright  Aeronautical 
Laboratories,  Wright- Patterson  AFB,  Ohio,  under  Project  2419,  "Non- 
metallic Structural  Materials",  Task  241903,  "Composite  Materials  and 
Mechanics  Technology",  Work  Unit  23190310,  "Durability  of  Composites 
and  Adhesives." 

The  work  reported  herein  was  performed  during  the  period  1 May  1976 
to  31  December  1977.  Dr.  H.  Thomas  Hahn  (AFML/MBM)  was  the  project 
engineer  and  Dr.  J.  Erikson  was  a visiting  scientist  from  the  National 
Defense  Research  Institute,  Stockholm,  Sweden,  for  the  period  15  June  1976 
to  31  August  1976. 

The  authors  wish  to  acknowledge  R.  Esterline  and  R.  Cornwell  of  the 
University  of  Dayton  Research  Institute  for  the  preparation  and  testing  of 
composite  tubular  specimens. 
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SECTION  I 


INTRODUCTION 

Many  theories  are  available  for  representation  of  the  strengths  of 
unidirectional  composites,  as  surveyed  in  [ 1 "J  and  [2],  As  regards  the 
matrix/interface- controlled  strength,  which  is  the  topic  of  discussion  in  the 
present  report,  these  theories  can  be  classified  by  the  degree  of  interaction 
between  the  normal  and  shear  stresses  and  also  between  tension  and  com- 
pression. For  example,  the  modified  von  Mises-Hill  criterion  [l]  takes 
into  account  the  coupling  between  the  transverse  normal  and  shear  stresses; 
however,  it  requires  a different  polynomial  when  the  transverse  stress  is 
compressive.  The  maximum  stress  or  maximum  strain  criterion  incorporates 
no  coupling  at  all;  it  is  based  on  the  independence  of  the  failure  modes.  The 
tensor  polynomial  criterion  [ 1 "J , however,  recognizes  full  coupling  and 
employs  only  one  polynomial. 

All  of  the  theories  have  their  advantages  and  disadvantages.  However, 
one  or  all  of  the  following  three  factors  can  be  cited  as  a reason  for  pre- 
ferring one  theory  to  the  others:  physical  foundation,  goodness  of  fit  with 
the  data,  and  convenience. 

In  most  cases,  the  goodness  of  fit  tests  have  been  performed  on  the 
off-axis  strength  without  much  success;  the  off-axis  strength  decreases 
rapidly  with  increasing  off-axis  angle,  thus  contributing  to  the  possiblity  of 
visual  deception  in  graphical  comparisons.  Consequently,  it  is  difficult  to 
distinguish  one  theory  from  another  graphically. 

Convenience  depends  on  the  type  of  application.  For  example,  the 
tensor  polynomial  criterion  is  simpler  in  digital  applications  whereas  the 
maximum  strain  criterion  is  more  convenient  in  graphical  applications. 

As  for  the  physical  foundations  behind  the  failure  theories,  one  may 
observe  that  the  von  Mises-Hill  criterion  assumes  failure  to  be  independent 
of  hydrostatic  stress  and  that  the  maximum  stress  or  maximum  strain 
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criterion  is  based  on  the  noninteraction  among  the  failure  modes.  However, 
no  physical  foundation  has  been  provided  for  the  tensor  polynomial  criterion. 

The  objective  of  the  present  report  is  to  discuss  some  of  the 
physical  evidences  in  support  of  the  tensor  polynomial  criterion  and  propose 
a corresponding  characterization  procedure  for  the  combined  loading  strength 
including  scatter.  The  discussion  is  limited  only  to  unidirectional  polymer 
matrix  composites  subjected  to  transverse  normal  and  longitudinal  shear 
stresses  so  that  composite  failure  is  controlled  by  the  matrix/interface 
properties.  Thus  the  appropriate  polynomial  is  of  the  form 


FZaZ  + FZZaZ  + F66<r6 
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The  foregoing  polynomial  follows  from  the  more  general  tensor  poly- 
nomial in  [ 1 ] in  the  absence  of  the  longitudinal  stress  in  the  fiber  direction, 
o^.  However,  recognizing  the  difference  between  the  fiber- controlled  fail- 
ure and  the  matrix/interface- controlled  failure,  we  propose  that  the  reduced 
polynomial  be  used  to  describe  the  matrix/interface- controlled  strength  even 
when  (Tj  / 0.  It  goes  without  saying  that  the  fiber- controlled  strength  is  then 
described  by 


F <r  + F cr 

11  111 
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SECTION  II 

PHYSICAL  BACKGROUND 


INTERRELATIONSHIP  BETWEEN  STRENGTH  AND  FRACTURE 
TOUGHNESS 


Just  like  homogeneous  brittle  materials,  unidirectional  composites 
exhibit,  when  the  artificially  introduced  crack  is  parallel  to  the  fibers  [3], 
those  fracture  characteristics  that  are  amenable  to  the  linear  elastic  frac- 
ture mechanics  predictions.  In  addition,  the  transverse  or  shear  strength 
of  unnotched  composite  is  known  to  depend  strongly  on  the  inherent  defects 
such  as  voids  and  interfacial  debonds  [4,5,6].  Thus,  if  these  defects  are 
regarded  as  typical  cracks,  the  unnotched  strength  can  be  predicted  from 
the  fracture  toughness  and  the  size  of  the  defect  [7,8]. 


Suppose  the  unnotched  strength  is  represented  by  the  polynomial, 

2 


2 

F c + F (t  + F cr 

2 2 22  2 66  6 
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(1) 


where  a is  the  transverse  normal  stress  and  cr  the  longitudinal  shear  stress, 

i.  6 

The  F's  are  the  components  of  strength  tensors. 


Now,  if  the  failure  is  assumed  to  initiate  at  an  inherent  crack  of  half 


length  a^  parallel  to  the  fibers.  Figure  1,  Equation  (1)  can  be  rewritten  in 


terms  of  the  nominal  stress  intensity  factors  k^ , k', 

d O 


A2k2  + A22‘'2 


k'Z  + A 


66k6 


(2) 


where 


k2 


<r  ,Va~ 

Z v o 


(3) 

(4) 


The  parameters  A's  are  related  to  the  F's  through 
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A2 

= F / V"a~  , 

Z o 

(5) 

A22 

= F /a  . 

ZZ  o 

(6) 

A66 

= F„/a 

66  o 

(7) 

In  the  foregoing  derivation,  the  surface  of  the  inherent  crack  is  not 

necessarily  normal  to  cr^,  Figure  1.  However,  when  is  tensile,  the 

inherent  crack  can  be  assumed  normal  to  and  the  nominal  stress  intensity 

factors,  k^,  k^,  reduce  to  the  mode  I and  mode  II  stress  intensity  factors, 

k , k^,  respectively.  Consequently,  Equation  (2)  can  be  considered  in  the 

first  quadrant  as  representing  a mixed-mode  fracture  toughness  envelope 

under  in-plane  loadings.  Moreover,  through  the  use  of  the  uniaxial  strengths 

and  corresponding  fracture  toughnesses  a can  be  determined  fi  ,o 

o 


or 


(8) 


a 

o 


(9) 


where  X^  and  X^  are  the  transverse  tensile  and  shear  strengths,  respectively, 
and  k and  k^  are  the  mode  I and  mode  II  fracture  toughnesses,  respectively. 

The  foregoing  hypothesis  can  be  checked  against  the  data  for  Scotchply 
1002.  The  appropriate  strength  and  fracture  properties  are  listed  in  Table 
1 [9].  The  inherent  crack  half  length  is  then 

a = 2.61  mm  or  3.53  mm  (10) 

depending  on  whether  Equation  (8)  or  (9)  is  used. 

There  are  many  factors  that  can  contribute  to  the  difference  in  the 
calculated  values  of  a^:  test  method,  material  variability,  and  the  assumption  of 
through-the- thickness  crack.  Although  effects  of  the  last  two  factors  cannot 


be  clearly  defined  at  present,  the  first  seems  to  point  in  the  right  direction. 

That  is,  k„  was  determined  from  a cantilever  beam  subjected  to  a concen- 
2c  J 

trated  load  [ 3]  . Since  this  test  is  closer  to  the  short  beam  shear  test  than 
to  the  off-axis  tension,  the  interlaminar  strength  should  be  used  for  X . 

Since  the  interlaminar  shear  strength  is  usually  higher  than  the  inplane  shear 
strength  that  was  used  in  Equation  (9),  the  resulting  a^  will  be  smaller  than 
3.53  mm  and  hence  closer  to  what  is  predicted  by  Equation  (8). 

For  the  foregoing  reasons  and  also  since  the  overall  comparison 
between  the  strength  and  fracture  toughness  under  combined  stress  is  of  inter- 
est, we  take  the  average  value 

a = 3.07  mm  (11) 

o 

and  proceed  to  investigate  the  consequences.  The  fracture  toughness  tensor 
components  follow  upon  substitution  of  the  F's  in  Table  1 and  a into  Equations 
(5  ) and  (7  ) : 


A2  = 

0.7735 

( MNm" 

3/2^" 1 

(12) 

A22  = 

0.1179 

(MNm~ 

3/2^ 

(13) 

A66  - 

0.07459 

(MNm 

3/2)"2 

(14) 

The  resulting  fracture  toughness  envelope  is  seen  to  agree  well,  in  the 
first  quadrant,  with  the  experimental  data  in  Figure  2.  It  should  be  noted 
that  the  data  were  obtained  from  the  artificially  introduced  cracks  which  are 
normal  to  cr^  [ 3]  . 

The  mode  II  fracture  toughness  increases  as  transverse  compression 
is  applied  because  of  the  friction  between  the  crack  surfaces  [9  ] . This 
increase  is  correctly,  although  lacking  in  quantitative  accuracy,  predicted 
by  the  polynomial  strength  criterion, but  not  by  the  other  criteria.  Also,  the 
coupling  between  kl,  and  k^  in  the  first  quadrant  can  be  described  neither  by 


the  maximum  stress  criterion  nor  by  the  maximum  strain  criterion. 

Thus,  it  can  be  concluded  that  the  polynomial  strength  criterion  agrees, 
when  the  transverse  stress  is  tensile,  with  the  basic  characteristics  of  the 
available  fracture  toughness  envelope.  As  the  transverse  compression  in- 
creases, however,  the  failure  will  in  general  start  from  an  inherent  crack 
whose  surfaces  are  inclined  to  the  applied  stress.  How  this  critical  angle  of 
inclination  depends  on  the  applied  stresses  should  be  known  if  the  entire 
failure  envelope  is  to  be  related  to  the  fracture  toughness.  By  way  of  illus- 
trating this  dependency,  we  discuss  the  energy  release  rate  approach  in  the 
following  subsection. 

2.  APPLICATION  OF  ENERGY  RELEASE  RATE  APPROACH 

Since  the  crack  extension  in  unidirectional  composites  is  at  least  macro- 
scopically  self- similar , the  energy  release  rate  can  be  calculated  easily. 

That  is,  for  a crack  parallel  to  the  fibers,  the  total  energy  release  rate 

1 

becomes  [10] 


G 

= G1 

+ G2  + g3 

(15) 

Here  the  energy  release  rates 

Gl'  G2 ’ 

G^  resulting  from  k^,  k2  and  k^. 

respectively,  are  given  by 

G1  = 

'Blkl  ’ 

B1 

1 

= Q eI 

T 

(16) 

/E  \1/2 

G2  = 

nB2kZ  ’ 

B2 

■ >.R) 

(17) 

G3  = 

nB3k3  ’ 

B3 

1 

1 / 2 

2(G  G _)A/ 

TT  LT' 

(18) 

Note  that  k^  is  the  mode  III  stress  intensity  factor.  The  orthotropy  correction 
factor  a is  defined  by 


I 


A 
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where  E is  the  longitudinal  modulus  parallel  to  the  fibers,  E the  transverse 
J -j  X 

modulus,  V the  minor  Poisson's  ratio,  G the  longitudinal  shear  modulus, 
1 L Xj  X 

and  the  transverse  shear  modulus. 

TT 

We  now  consider  a crack  of  length  2a  which  is  not  normal  to  but  in- 

o 2 

dined  by  an  angle  G,  Figure  1.  In  terms  of  the  applied  stresses  cr  and  a , 

L*  6 

the  stress  intensity  factors  for  this  crack  are  expressed  as  follows; 


kl 

(. 

- o-.va  cos  0 , a_  A o 

2 o 2 

(20) 

kl 

= 0 , cr^  < 0 

(21) 

k2 

= <r.  v/a  cos  0 

6 o 

(22) 

k3 

= sin  0 cos  0 

2 v o 

(23) 

Note  that  in  the  calculation  of  k^  no  friction  has  been  taken  into  account. 

When  cr^  is  positive,  G becomes  stationary,  i.e.,  dG/dG  = 0,  at  the 
following  angles: 


G = 0 or  rr/2 

. - 1 r i bi  + b2(°v°2) 

0 = sm  [z  - -(b3  - B“ 


.2  1 
a2] 

1>  . 


On  the  other  hand,  if  <r  is  negative,  then  G = 0 and  the  angle  given  by 
Equation  (25)  is  changed  to 

T 211/2 

B„  la,  \ 


_ . -1  1 “2  /a6\ 
G = sin  — - r— — [ — ) 

2 2B3\°z) 


Elastic  properties  of  commonly  used  composites  are  listed  in  Table  2 
together  with  the  calculated  values  of  B's.  For  .hose  composites,  we  have 
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(27) 


2B1>B3  (27) 

and  therefore  Equation  (25)  cannot  be  satisfied.  Thus,  when  is  tensile, 

the  maximum  energy  release  rate  occurs  at  9 = 0 independently  of  the  shear 

stress  a and 
6 

G = G | = 77a  (B  cr  + B_c r ) (28) 

max  1^  _ o 1 2 2 6 

However,  when  is  compressive,  the  critical  crack  orientation  depends 

on  the  stress  ratio  <r  /cr  . That  is,  if  Icr  /<r  I > (B_/B_)^  then 

o £ 1 3 £ 


G = G 

max 


77a  B_cr. 
o 2 6 


1 /2 

On  the  other  hand,  if  cr  / cr  < (B  /B  ) , then 

6 c 3 L 


G = G 

max 


2 1 J_  ^3/^2  \ _1  ^2 /^6  \ 

o 2a6  2 + 4 B2\cr6  ) + 4 B3(<r2/ 


where  9q  satisfies  Equation  (26).  Note  that,  when  <7^=0,  i.e.,  pure  compression, 


one  obtains 


G = rra  B_<r_/4 

max  o 3 2 


If  the  energy  release  rate  is  applicable,  the  uniaxial  strengths  are 
related  to  the  critical  energy  release  rate  G^  by 


77a  B_ 
o 2 


• M • 


fz  /M 

x6  " v bi  y 


/M 

W 


1/2 


1/2 


= 21 


(34) 


(35) 


where  X^  is  the  compressive  transverse  strength.  The  combined- stress  fail- 
ure criteria  then  follows  from  Equations  (28)-(30).  That  is,  when  cr^  > 0, 
Equation  (28)  reduces  to 


(i)‘  •(£)' 


= 1 


(36) 


When  a < 0,  on  the  other  hand,  two  different  criteria  follow  depending  on  the 

^ 1/2 

ratio  Jo-  /a  |.  If  | cr  / cr  | > (B  /B_)  , then  the  shear  strength  is  independ- 


ent of  c 


V 


2' 

1 • 6 • 


= x. 


(37) 


.1/2 


However,  if  | cr  /cr  | < ( B /B_)  , then  the  failure  criterion  is  given  by 


= l 


(38) 


Table  3 lists  the  appropriate  strength  ratios  for  various  composites, 
both  experimental  and  predicted.  Insofar  as  the  ratio  X^/X^  is  concerned, 
the  energy  release  rate  approach  seems  to  yield  a fairly  good  correlation  with 
the  data  for  both  graphite/epoxy  composites.  However,  the  correlation  is 
very  poor  for  the  ratio  X^/X^  irrespectively  of  the  type  of  material. 

In  any  case,  the  energy  release  rate  implies  a strong  interaction 

between  <j  and  cr  . Aside  from  the  applicability  of  the  energy  release  rate 
2 6 

approach,  there  are  many  assumptions,  such  as  the  noninteracting  through- 
the-thickness  cracks,  the  mode  III  fracture  in  compression,  etc.,  that  have 


- 


to  be  validated  in  order  to  account  for  the  discrepancy  between  the  theory  and 
the  data.  Still,  what  is  interesting  is  that  the  dependence  on  the  inherent 
crack  length  of  the  strength  under  combined  stresses  simply  follows  from  the 
failure  envelope  if  cr^  and  u ^ are  replaced  by  <r2Y^a~  anc*  a£,V  aQ  > respectively. 
Thus,  it  is  plausible  to  assume  that  Equations  (1)  through  (7)  are  valid  in  the 
entire  sPace  anc^  that  the  A's  are  independent  of  a^. 

3.  DERIVATION  OF  STRENGTH  DISTRIBUTION  FROM  CRACK  LENGTH 
DISTRIBUTION 


Now  that  the  relationship  between  the  strength  and  the  inherent  crack 
length  is  known,  we  can  derive  a strength  distribution  under  combined  loading. 
To  this  end  we  further  note  that  the  critical  crack  orientation  is  independent 
of  the  crack  length,  as  shown  in  the  preceding  subsection.  Thus,  we  can 
assume  that  the  scatter  in  strength  is  solely  due  to  the  variation  in  the  crack 
length. 


Suppose  the  crack  half  length  has  the  cumulative  distribution 


P(a) 


(39) 


where  a and  o are  the  characterization  parameters.  Since  the  A's  in  Equations 
(5)  - (7)  can  be  considered  as  deterministic  material  constants,  we  obtain  the 
following  distributions  of  F's: 


P(F2)  = exp 


p(F22)  = exp 


p(F66)  = exp 


■“  /£  \q/2-i 

(ft) 

r /F  \a/Zn 

(ft) 


A , A 

where  F's  correspond  to  a through  Equations  (5)  - (7). 


(40) 


(41) 


(42) 
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We  now  introduce  the  concept  of  a strength  vector  [9]  . Suppose  failure 
occurs  at  the  stresses  and  a^>  The  pair  (a^,  a^)  then  defines,  in  the  stress 
space,  a vector  which  is  called  the  strength  vector.  Since  the  failure  envelope 
represents  a set  of  pairs  of  failure  stresses,  any  point  on  the  failure  envelope 
has  a strength  vector  associated  with  it.  From  Equations  (1)  and  (5)-(7),  the 
magnitude  a of  this  strength  vector  (<r2,<r6),  called  the  combined  strength  here- 
after, is  shown  to  be  proportional  to  1 /%/  a : 


2 2 2 1/2 
[(A^  + 4A22)cos  (t>  + 4A6fesin  0 ] - A2cos 

~ ' 2 ~2 

2(A22cos  0 + A^sin  <M 


where 


0 = tan  (cr^/(r2^ 


Then,  introducing  the  characteristic  combined  strength, 

A i 

a = o- 

I A 


we  can  rewrite  Equation  (43)  as 


,A  .1/2 

(!) 


Therefore,  the  probability  of  the  combined  strength  being  greater  than  <r, 
R(cr),  is  expressed  in  the  following  simple  form: 


R(cr)  = exp 


[-(*)“] 


Note  that  the  characteristic  combined  strength  a varies  with  the  angle 
0.  However,  the  distribution  of  the  ratio  a/a  remains  independent  of  0. 
Also,  a/a  has  the  same  distribution  as  do  the  uniaxial  strengths.  In  fact. 
Equation  (47)  reduces  to 
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As  is  clear  by  now,  the  shape  parameter  Q is  independent  of  the  mode 
of  loading.  Table  4 lists  values  of  a for  some  composites  available  in  the 
literature.  No  definitive  relationship  is  seen  between  a and  the  mode  of  load- 
ing although  some  of  the  inconsistencies  may  be  due  to  the  difference  in  the 
test  methods  employed.  Therefore,  in  the  following  section,  we  assume  a is 
independent  of  the  type  of  loading  and  adopt  the  procedure  based  on  tr  to  char- 
acterize the  combined  strength  data. 


i 
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SECTION  III 

ANALYSIS  OF  COMBINED  LOADING  DATA 


1.  CHARACTERIZATION  OF  STRENGTH 

Strength  data  were  obtained  by  testing  off-axis  tubes  in  combined  axial 
and  torsional  loadings.  The  experimental  procedures  are  described  in  f 17]. 

The  strengths  listed  in  Table  5 are  fit  by  the  equation 


Fq2°2  + Fo22°2  + Fo66°6 


(51) 


where  F^'s  are  determined  by  the  least  squares  method.  That  is.  Equation 

(51)  for  each  set  of  data  (cr^,  <r^),  is  rearranged  as  follows 

2 b 


> 

Fo2 

1 

42) 

42,z 

< 

Fo22 

1 = < 

1 

- 

- 

- 

Fo66 

L J 

- 

i 

q 

4n’2 

^ 1 > 

(52) 


Using  an  abridged  notation  for  Equation  (52),  i.e. 

0<O  - (>)■ 


we  can  determine  F 's  from  the  following  equation: 
o 

-lr  -.T / 


(53) 


(f  ) = (OfM  r'H  (') 


(54) 


where  the  subscripts  T and  -1  stand  for  the  transpose  and  inversion,  respec- 
tively. The  results  are 


F _ = 3.  376  x 10"2  (MPa)”  1 

o2 


1 3 


(55) 


o22 

= 4.721  x 10'4  (MPa)"2 

(56) 

-4  - 2 

066 

= 2.  384  x 10  (MPa) 

(57) 

Next,  in  order  to  determine  F's,  a nondimens ional  strength  parameter 
s is  defined  by 


s = a/cr 

o 


(58) 


where  a is  the  actual  combined  strength  and  a ^ is  the  combined  strength  pre- 
dicted by  the  failure  envelope.  Equation  (51).  That  is,  if  (a^.o^)  is  a pair  of 
stress  components  at  failure,  then 


/ 2 Z\1'2 

* = (a2  +<r6) 


(59) 


(ao2+(yo^ 


2 M/2 


where  (cr  , cr  .)  satisfies  Equation  (51)  and 
o2  oo 


(60) 


= s 


(61) 


o2  06 

Recalling  that  <r  has  a Weibull  distribution,  Equation  (47),  we  deduce 


R(s)  = exp  -I- 


[•(-;)“] 


A 

A cr 

s = — 
cr 


(62) 


The  experimental  values  of  s are  plotted  in  Figure  3 according  to  the 
linearized  format  of  Equation  (62); 


A 


f n(  - t nR)  = q f n s - a!n  s 
The  median  rank  was  used  for  R; 


R = 1 


0.3 


N + 0.4 


(63) 


(64) 
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where  j is  the  ordinal  number  of  strength  and  N the  total  number  of  data. 

The  data  represented  by  the  closed  symbols  are  seen  to  deviate  much 
from  the  trend  exhibited  by  the  open  symbols.  Two  of  those  four  data  were 
obtained  from  badly  misaligned  15-degree  off-axis  tubes;  therefore,  these  can 
be  regarded  as  invalid  data.  The  remaining  two  did  not  exhibit  any  apparent 
anomalies.  However,  it  is  still  possible  that  a slight  misalignment  may  have 
caused  premature  failure  since  these  two  are  90-degree  tubes. 

Therefore,  the  aforementioned  four  points  are  discarded  and  new  F s 

o 

are  calculated  with  the  results 


Fo2  = 

3. 126  x 10~2 

(MPa)'1 

(65) 

„ -4 

-2 

Fo22 

4.428  x 10 

(MPa) 

(66) 

„ „ -4 

-2 

Fo66 

2. 391  x 10 

(MPa) 

(67) 

The  distribution  of  s based  on  the 

newly  determined  F ' s 

o 

is  shown  in 

Figure 

The  best-fit  curve  results  from 


a = 4.745  , s = 0.9695  (68) 


and  the  correlation  coefficient  is  0.9922.  Furthermore,  the  corresponding 

A 

F's  are  obtained  from  the  following  equations: 


(69) 


a 2 
s 


(70) 


Thus,  the  matrix/interface- controlled  strength  of  the  composite  studied 

A 

is  characterized  by  the  distribution  of  < r,  Equation  (47)#  and  by  the  F's: 


A 

F2 

= 

3.224  x 

IQ'2 

(MPa)-1 

(71) 

A 

F22 

= 

4 . 567  x 

1C’4 

(MPa)'2 

(72) 

A 

F66 

= 

2 . 466  x 

io'4 

(MPa)'2 

(73) 

Now,  given  any  applied  stresses,  the  survival  probability  can  be  calculated. 

As  an  example,  suppose  the  following  stresses  are  applied  to  the  pre- 
viously characterized  composite: 

cr2  = -40  MPa  , v = 60  MPa  (74) 

The  stress  ratio  cr/cr  follows  from  Equations  (71)-(73)  as 

| = °-779  (75) 


Thus,  the  probability  of  survival,  R,  is  given  by 

R = exp  [-(0.779)4'745]  = 0.737  (76) 

2.  SIZE  EFFECT 

To  see  if  the  weakest  link  theory  based  on  s can  account  for  the  size  effect, 
some  data  were  obtained  from  90- degree  tubes  having  gage  length  reduced  to  3/8 
the  original  gage  length.  Described  in  the  following  are  the  results,  both 
analytical  and  experimental. 


Following  the  same  procedure  as,  e.g.,  in  [18],  but  using  s,  we  obtain 
the  ratio  of  the  characteristic  strength  parameters  as 


A 

s 


(!) 


1/4.745 


1.2296 


(77) 


Here  the  primed  quantities  are  associated  with  the  short  tubes  and  V is  the 
, A 

volume.  The  F's  for  the  short  tubes  are  then 
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i 


A 

F ' 

2' 

= 

6 A,  A, 

F2  s/s'  = 

2 . 542  x 

io'2 

(MPa)"  1 

(78) 

A 

F 1 
22 

= 

2.9  39  x 

io“4 

-2 

(MPa) 

(79) 

A 

F 1 
66 

= 

^6^’  = 

1.581  x 

io'4 

-2 

(MPa) 

(80) 

The  resulting  failure  envelope  is  compared  with  the  experimental  results 
in  Figure  5.  The  broken  curve  is  the  best  fit  of  the  data.  Although  any  defini- 
tive conclusions  require  more  data,  the  comparison  seems  quite  encouraging. 


SECTION  IV 
CONCLUSIONS 

We  have  presented  a method  of  characterizing  the  matrix/interface- 
controlled  strength,  including  scatter,  of  unidirectional  composites  under 
combined  transverse  normal  and  shear  loading.  It  was  assumed  that  failure 
initiates  at  inherent  cracks  parallel  to  the  fibers  and  that  the  scatter  in 
strength  is  a manifestation  of  nonuniform  crack  length  having  a certain  dis- 
tribution. The  energy  release  rate  approach  was  discussed  as  a guide  for  a 


possible  relationship  between  the  combined  strength  and  the  inherent  crack 
length. 

The  proposed  failure  criterion  is  a second  order  polynomial.  The 
scatter  is  described  by  the  strength  vector  whose  magnitude  has  a Weibull 
distribution.  The  method  holds  promise  as  a means  of  accounting  for  the  size 
effect  under  combined  loading. 

In  the  energy  release  rate  approach,  actual  defects  were  replaced  by 
the  noninteracting  through-the-thickness  cracks  and  the  critical  crack  orien- 
tation was  assumed  to  depend  on  the  applied  stresses  in  a deterministic 
manner.  The  relaxation  of  the  first  assumption  calls  for  a three-dimensional 
stress  analysis  and  the  second  assumption  can  be  alleviated  by  employing  a 
statistical  approach,  as  was  done  for  homogeneous  materials  in  [19].  How- 
ever, a major  obstacle  to  the  establishment  of  a strength- fracture  toughness 
relationship  seems  to  be  the  lack  of  an  appropriate  fracture  criterion  under 
combined  state  of  stresses. 
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TABLE  1 


STRENGTH  AND  FRACTURE  PROPERTIES 
OF  SCOTCHPLY  1002(a) 


X2(MPa) 

X (MPa) 

D 

- 3/2 

k.  (MNm  ) 

1 c ' 

- 3/2 

k2c(MNm  ) 

20.0 

66.2 

1.022 

3.934 

F £ (GPa)'1 

F22  (GPa)-2 

F66  (GPa,-2 

42.79 

360.76 

228.24 

(a)  Data  taken  from  [9] 
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TABLE  2 


ELASTIC  PROPERTIES 

Gl/Ep 

Scotchply  1002 

Gr/Ep 

T300/5208 

Gr/Ep 

AS/3501 

B/Ep 

Narmco  5505 

EL(GPa) 

34.47 

181. 33 

137.69 

207.53 

ET(GPa) 

11.49 

10. 34 

9.65 

18.82 

GLT(GPa) 

4.86 

7.17 

4.21 

5.24 

VLT 

0.05 

0.28 

0.  30 

0.21 

GTT(a)(GPa) 

4.  32 

3.89 

3.63 

7.08 

BjlMPa)’  1 

81.24 

66.43 

86.  35 

54.  16 

B2(MPa)_1 

46.90 

15.86 

22.86 

16.  31 

B3<  MPa)’  1 

109.12 

94.69 

127.88 

82.11 

Reference 

[3] 

[11] 

[12] 

[13] 

(a)  Gtt  = Et/[2(1  + vTT)]  , vTT  = 0.33. 


STRENGTH  RATIOS 
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STRESS  COMPONENTS  AT  FAILURE 


(a)  Retested  after  failure. 
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Figure  3.  Linearized  Plot  of  s. 
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TRANS.  S 


Figure  5. 


Predicted  Failure  E 
for  Short  Tubes. 
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